a-nucleus potentials, a-decay half-lives, and shell closures for superheavy nuclei 
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Systematic a-nucleus folding potentials are used to analyze a-decay half-lives of superheavy nuclei. 
Preformation factors of about several per cent are found for all nuclei under study. The systematic 
behavior of the preformation factors and the volume integrals of the potentials allows to predict 
a-decay energies and half-lives for unknown nuclei. Shell closures can be determined from measured 
a-decay energies using the discontinuity of the volume integral at shell closures. For the first 
time a double shell closure is predicted for Zmaglc = 132, -/V mag ic = 194, and j4 mag ic = 326 from the 
systematics of folding potentials. The calculated a-decay half-lives remain far below one nanosecond 
for superheavy nuclei with double shell closure and masses above A > 300 independent of the precise 
knowledge of the magic proton and neutron numbers. 
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The a-decay of superheavy nuclei has been studied in- 
tensively in the last years 0, 0, 0, II 0, 0, 0, 0, ©, G3 • 
In many papers a simple two-body model was applied 
and in most papers a potential was derived which 
was able to fit the measured a-decay half-lives of the an- 
alyzed nuclei._ However, most of the studies (with the 
exception of 0) did not attempt to use these potentials 
for the description of other experimental quantities like 
e.g. a scattering cross sections or (n,a) or fusion reaction 
cross sections. 

Therefore, an alternative approach was followed in [l2l ]. 
Now the simple two-body model has been combined with 
systematic a-nucleus folding potentials which are able to 
describe various properties, and the ratio between the 
calculated half-life Ty}f a and the experimental half-life 
T°J^ a has been interpreted as preformation factor P of 
the a particle in the decaying nucleus. Besides a sys- 
tematic behavior of the volume integrals of the folding 
potentials, preformation factors of the order of a few per 
cent were found for a large number of nuclei [H, EH • 

Only for very few light nuclei some levels have been 
found where a simple two-body model can exactly repro- 
duce the experimental half-lives or widths, e.g. for 6 Li 
= 2 H ® a 114] or for 8 Bc = a <8> a 0, Already 
for 20 Ne = 16 <g> a the calculated widths are somewhat 
larger than the experimentally observed ones 16]. Any 
simple two-body model with a realistic potential must 
provide half- lives identical or shorter than the experimen- 
tal value, because the two-body model assumes a pure a 
cluster wave function by definition, whereas any realis- 
tic wave function is given by the sum over many different 
configurations. Thus, preformations of a few per cent are 
a quite natural finding for superheavy nuclei. 

The following ingredients were used in this paper. The 
a-nucleus potential was calculated from a double-folding 
procedure with an effective interaction [l2|, [I3> EH ■ The 
nuclear densities were taken from for the a particle 
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with properly scaled radius 



parameter r ~ A T . The total potential is given by 
the sum of the nuclear potential Vjy(r) and the Coulomb 
potential Vc{r): 

V(r) = V N (r) + V c (r) = A V>(r) + V c (r) (1) 

The Coulomb potential is taken in the usual form of a 
homogeneously charged sphere where the Coulomb ra- 
dius Rc has been chosen identically with the rms radius 
of the folding potential Vf , and the folding potential Vf 
is scaled by a strength parameter A which is of the or- 
der of 1.0 — 1.3. This leads to volume integrals of about 
Jfj w 300MeVfm 3 for all nuclei under study and is in 
agreement with systematic a-nucleus potentials derived 
from elastic scattering [U HH, H [H, [H [H [H . (Note 
that as usual the negative sign of Jr is omitted in this 
work.) Bound state properties of 212 Po = 208 Pb ® a have 
been analyzed successfully using the same potential [13] ■ 
The centrifugal potential has been omitted for L = de- 
cays. The following study is restricted to even-even nuclei 
because the additional centrifugal barrier may influence 
the a-decay half-life for decays with L ^ 0. 

The quotations of the volume integral Jr and the po- 
tential strength parameter A are practically equivalent for 
this paper. If one wants to compare this work to folding 
potentials with a different nucleon-nucleon interaction or 
even to potentials with a different parametrization (e.g., 
Woods-Saxon potentials), the volume integral Jr has to 
be used. Therefore the following discussion is restricted 
to volume integrals. Nevertheless, most figures provide 
both quantities Jr and A. 

The a-decay width T Q is given by the following formu- 
lae [U: 



PF — exp 
4 M 



rr 3 

-2 / k(r)dr 

J r 2 



(2) 



with the preformation factor P, the normalization factor 
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and the wave number k(r) 
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[x is the reduced mass and E is the decay energy of the 
a-decay which was taken from the computer files based 
on the mass table of Ref. [28| or from Table 1 of Q. 
The Ti are the classical turning points. For + — > + 
s-wave decay the inner turning point is at n = 0. ri 
varies around 9 fm, and r% varies strongly depending on 
the energy. The decay width T a is related to the half-life 
by the well-known relation T a = 7im2/Tx/2 )Q! . Following 
Eq. ([2]), the preformation factor may also be obtained as 
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where T, 



or TffiV are calculated from Eq. © with P 



1/2, a 



1. For completeness, I define the here predicted half-life 
for unknown nuclei as T, p !° —T?fl c IP. 

1/2. a 1/2, a/ 

The potential strength parameter A was adjusted to 
the energy of the a particle in the a emitter (A + 4) = 
A (g> a. The number of nodes of the bound state wave 
function was taken from the Wildermuth condition 



= 2N + L = J2( 2n i + l i)=^2> 



(6) 



where Q is the number of oscillator quanta, N is the 
number of nodes and L the relative angular momentum 
of the a-core wave function, and qi — 2rii + li are the 
corresponding quantum numbers of the nucleons in the 
a cluster. I have taken q = 5 for 82 < Z, N < 126 
and q = 6 for N > 126 where Z and N are the proton 
and neutron number of the daughter nucleus. The above 
definition of Q slightly deviates from the semi-classical 
Bohr- Sommerf eld quantum number G as mostly used. 
One finds G w 22.5 for all nuclei with Q = 22. 

Various attempts have been made to determine the 
preformation factors P experimentally or theoretically 
[29l [30l l3ll ]. The usage of a simple two-body wave 
function in connection with the Wildermuth condition 
Eq. ^ is obviously a quite simple approximation for 
the description of the complex many-body wave func- 
tion of a superheavy nucleus which was analyzed e.g. in 
[13, HE HE HE Hil • Nevertheless, the preformation factor 
defined as ratio P = T^f 2 c jT^ a in Eq. © may be un- 
derstood as effective preformation factor. The obtained 
values for P do only show small variations and can thus 
be used for the prediction of half-lives of unknown su- 
perheavy nuclei in a consistent way. A full discussion of 
preformation factors is beyond the scope of the present 
paper. 

The resulting preformation factors P for even-even nu- 
clei are shown in Fig. [1] An average value of P « 8 % is 
found. Almost all results lie within a bar of uncertainty 
of a factor of three. This uncertainty is identical to the 
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FIG. 1: The preformation factors P are shown for several 
superheavy a emitters. The horizontal lines indicate an aver- 
age value of P w 8 % (full line) and typical uncertainties of a 
factor of three (dotted lines). 



results of 0, However, the values for P are much 
smaller in this work (see discussion above). A table of 
the results will be given in a forthcoming paper. 

There are two different ways in this simple two-body 
model to obtain larger a-decay half-lives Ty^ c a and thus 
larger preformation factors P as derived from Eq. l[5]). 
First, very narrow potentials can be used. In this case 
the attractive nuclear potential becomes negligible in the 
region of the Coulomb barrier thus effectively increas- 
ing the barrier and increasing the a-decay half-life. This 
idea was followed e.g. in and the differences to the 
systematic folding potential in the present work are il- 
lustrated in Fig. 1 of fl2| . A very narrow potential as 
used in [l(| is probably not able to describe quantities 
beside the a-decay half-life. Second, a smaller quantum 
number (G or Q) can be used. This idea was followed 
in [l|. In this case the attractive nuclear potential is re- 
duced at all radii, thus again effectively increasing the 
Coulomb barrier and the a-decay half-life. Many quan- 
tities are mainly sensitive to the tail of the wave func- 
tions at large radii outside the nuclear potential which 
leads to discrete ambiguities for the volume integral Jr 
of a- nucleus potentials (the so-called "family problem" ) . 
However, it has been found in the last years that system- 
atic a-nucleus folding potentials have volume integrals 
J R around 300 MeV fm 3 H3, E3, HI HI IM HE HI com- 
patible with the quantum number Q used in the present 
work and incompatible with the smaller G used in [l[. 

In principle, the application of a semi-classical model 
is not necessary for the calculation of a-decay half-lives 
or widths. From the potential in Eq. ([T]) one can directly 
calculate the wave function and the width of the decaying 
state. However, in practice this is difficult because of the 
low energies and extremely small widths of the states. 
For 212 Po = 208 Pb ® a such a fully quantum-mechanical 
calculation is possible at the limits of numerical stability. 
Fig. [2] shows the scattering phase shift Sl for the L = 
partial wave as a function of energy which is given by 
E = E Q + i x AE with E = 8.954088523002 MeV and 
AE = 2 x 10" 15 MeV. The points are the results of a 
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FIG. 2: Phase shift Sl for the L = partial wave for the 
system 212 Po = 208 Pb ® a. The derived width from Eq. © 
is T = 3.52 x 10~ 14 MeV. Note the extremely small stepsize 
of the calculation of AE = 2.0 x 10~ 15 MeV! 



phase shift calculation, the line is a fit to the points using 
the formula for narrow resonances 



Sl(E) = arctan 



r 



2(E R - E) 



(7) 



with E R = 8.9541 MeV and T = 3.52 x l(T 8 eV which 
translates to Tffi a = 13.0 ns. The semi-classical ap- 
proximation yields Ty^ a = 8.7ns which is about 30% 
lower than the precise quantum-mechanical value. A sim- 
ilar result is obtained for the decay of 8 Be where one 
finds r Q = 6.7 eV for the quantum-mechanical calcula- 
tion and r a = 10.5 eV for the semi-classical approxima- 
tion. As pointed out above, the preformation factors are 
close to unity for 8 Be with P = 100 % (65 %) for the 
quantum-mechanical (semi-classical) calculation and of 
the order of a few per cent for 212 Po with P = 4.3 % 
(2.9%). These results confirm the applicability of the 
semi-classical model within uncertainties of about 30%. 

It is interesting to use the systematic folding potentials 
for the prediction of properties of unknown superheavy 
nuclei like a-decay energies, a-decay half-lives, and shell 
closures above N,Z = 126. The basic building block is 
the smooth behavior of the strength parameter A of the 
folding potential and the resulting volume integrals Jr 
(see Fig. |3] and Table I of (H). 

Within one major shell, one finds variations of Jr from 
about Jr « 335 MeV fm 3 at the lower end of a shell to 
about Jr k 280 MeV fm 3 at the upper end of a shell. 
Between neighboring nuclei the variation in Jr is be- 
low AJr < 5MeVfm 3 . This allows first the determi- 
nation of a-decay energies for up to now unknown nu- 
clei. As an example, one finds for the decay of 298 120 
-> 294 118 a volume integral of Jr rj 296 MeV fm 3 cor- 
responding A = 1.138. This leads to a decay energy of 
E = 12.87 MeV. The a-decay half-life can be estimated 
using the given energy and an average preformation fac- 
tor of P « 8% leading to ^""w Sfis. Whereas the 
uncertainties for the volume integral Jr and the derived 
a-decay energy are small, the uncertainty of the a-decay 
half-life is strong because of the exponential energy de- 
pendence. For a potential strength enhanced (reduced) 
by 2% one finds the a-decay energy E = 10.70 MeV 
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FIG. 3: Volume integrals Jr for superheavy nuclei as a func- 
tion of Zd (upper), Nd (middle), and Ad (lower). Within a 
major shell one finds a smooth decrease of Jr with a minimum 
value around Jr ~ 280 MeV fm 3 . From the fits to the data 
points one can directly see that Jr drops below 280 MeV fm 3 
at2a 130, N m 192, and A w 322 leading to the magic 
numbers Z magic = 132, iV ma gic = 194, and A maglo = 326. The 
open circles are extrapolations for 294 118 and the lower limit 
of Jr > 279 MeV fm 3 close to the next closed shells (see text). 
The open squares are nuclei with Q = 18 ( 208 Pb = 204 Hg ® 

206 Pb ® a). 
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(E = 14.98 MeV) and T*% a = 0.97 s (T*% a = 1.6 ns) 
again using P = 8 %. A variation of the potential 
strength of 1 % corresponds to a variation of the a-decay 
energy of about 1 MeV which is comparable to the uncer- 
tainties of mass formulae [37j • As usual, the reliability 
of such an extrapolation decreases for nuclei with masses 
far above the heaviest known nuclei. However, the un- 
certainties for closed-shell nuclei remain small because of 
the well-defined volume integral Jr for such nuclei which 
can be studied at the shelle closures at N = 82, Z = 82, 
and N = 126. 

Shell closures can be seen as discontinuities in the vol- 
ume integrals, see Figs. |3] and [U Whereas the varia- 
tion between neighboring nuclei remains below AJr < 
5 MeV fm 3 , at shell closures one finds a strong increase of 
Jr which is directly related to the increase of the quan- 
tum number Q. Because shell closures are not known 
a priori for superheavy nuclei, Fig. [4] analyzes the vol- 
ume integrals around the shell closure at N — 82 for 
Xe, Ba, Ce and Nd isotopes. Below N — 82, the 
wave functions are characterized by Q = 16 (full black 
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FIG. 4: Potential strength parameter A (upper part) and vol- 
ume integrals Jr (lower) around the shell closure N — 82 for 
130 " 136 Xe ® a (diamonds), 132 ~ 142 Ba® a (circles), 134 ~ 144 Ce 
<g> a (triangles), and 140_146 Nd <g> a (squares) isotopes. One 
can clearly see the discontinuity at the shell closure N = 82 
(see text). 



symbols), and the volume integrals are slightly above 
J R = 280MeVfm 3 . Above N = 82 one finds volume inte- 
grals around 310 MeV fm 3 for wave functions with Q — 18 
(open symbols; see also [l2j )■ The small grey symbols are 
calculated above the shell closure at TV = 82 without in- 
crease of the quantum number Q. Here one finds low 
volume integrals significantly below Jr = 280 MeV fm 3 
which differ by more than AJr = lOMeVfm 3 from 
the neighboring values. The behavior of the potential 
strength parameter A is similar to Jr (see Fig.[4j. 

A similar behavior is found at the shell closures at 
Z — 82 and N = 126 around 208 Pb. Compared to neigh- 
boring values, the volume integrals Jr are reduced by 
more than A Jr = 10 MeV fm 3 if one neglects the increase 
of the quantum number Q, and the absolute values of the 
volume integrals drop below Jr = 280 MeV fm 3 . Conse- 
quently, changes in Jr by about lOMeVfm 3 or values of 
Jr significantly below 280 MeV fm 3 are clear indications 
for the crossing of a major shell. Because the determina- 
tion of the volume integral Jr requires only the knowl- 
edge of the a-decay energy, the measurement of one sin- 
gle quantity may be sufficient for the determination of a 
double closure: as soon as Jr drops below 280 MeV fm 3 , 
magic neutron or proton numbers have been crossed! 

This systematic behavior of a-nucleus potentials allows 
further a prediction of magic numbers in a yet unknown 
mass region above A > 300. The smooth energy depen- 
dence of the volume integrals Jr in Fig. [3] is fitted using 
a second-order polynomial for all nuclei with Q — 22 (full 
lines in Fig.©. Jr drops below 280 MeV fm 3 at Z = 130, 
N = 192, and A = 322 which means that the nucleus 
326 132 = 322 130 a is the heaviest nucleus which can be 
described using a potential with Jr > 280 MeV fm 3 and 
Q = 22. Increasing Z or N leads to Jr below its lower 



limit, and thus the magic numbers Z n 



132 ± 4, 



Vn 



194 ± 4, and A n 



326 ± 6 can be derived 



from Fig. [3] 

The a-decay half-life of the doubly-magic nucleus with 
Zrnagic = 132, 7V magic = 194, and A magic = 326 can be 
calculated using the volume integral Jr = 279.2 MeV fm 3 
(taken from 208 Pb = 204 Hg ® a). One finds the energy 
E = 18.26 MeV and the corresponding half-life T^f 2 c a 
— 1.16 x 10 -12 s with P = 1. Again using P = 8 %, a re- 
alistic prediction of the half-life is a = 1-5 x 10~ n s. 
Including the uncertainty of P, the half-life remains be- 
low 10 -10 s. The uncertainty of the volume integral Jr 
at closed shells is smaller than 1 %. Increasing Jr by 
1 % reduces the a-decay energy by about 1 MeV and in- 
creases the a-decay half-life by about a factor of 20. In 
any case, the half-life remains below 1 ns. 

The lower limit of Jr has also been applied for the 
calculation of the half-life of 310 126 -> 306 124 with the 
widely discussed shell closures at Z = 126 and N = 184 
(e.g., [38|,l39|); but also other magic numbers have been 
discussed (e.g. [4(| )• Here one obtains the a-decay energy 
E = 18.82 MeV and the corresponding a-decay half-life 
T i/2° a = 2- 1 x 10~ 14 s (P = 1). The realistic prediction 
using P = 8 % is a = 2.6 x 10~ 13 s, and including all 

uncertainties the half-life remains far below 10~ n s. The 
above calculations indicate clearly that one cannot expect 
that any superheavy nucleus above A > 300 with magic 
proton and neutron numbers (whatever these numbers 
are) has a half-life significantly above one nanosecond. 

Because of the significant variation of the decay en- 
ergy E from about 4 MeV to about 12 MeV for known 
superheavy nuclei and up to about 20 MeV for the pre- 
dicted but yet unknown doubly-magic superheavy nu- 
clei one may expect a correlation between the potential 
strength parameter A or the volume integral Jr and the 
decay energy E. This relation is analyzed in Fig. [5] for 
superheavy nuclei and Fig. [5] for nuclei around N = 82. 

Fig. [5] seems to indicate that larger decay energies E 
are correlated to smaller volume integrals Jr. However, 
the underlying reason for this energy dependence of Jr 
is the smooth variation of Jr within a major shell (see 
above). At very small energies one finds again a small 
volume integral of J R « 280 MeV fm 3 for 208 Pb = 204 Hg 
(g) a. As can also be seen from Fig. [HI the volume integrals 
do not depend on the energy E: above TV = 82 one finds 
J R w 310 MeV fm 3 for bound (E < 0) and unbound (E > 
0) nuclei, and below N = 82 one finds J R » 280 MeV fm 3 , 
again for bound and unbound nuclei. 

In conclusion, systematic folding potentials can be 
used for the calculation of a-decay half-lives of super- 
heavy nuclei. Additionally, the systematic behavior of 
the volume integrals allows to predict a-decay energies 
and half-lives of yet unknown nuclei. The magic num- 
bers Z magic = 132, iVmagic = 194, and A magic = 326 have 
been derived from the discontinuities of the volume in- 
tegrals at shell closures. There is strong evidence that 
a-decay half-lives remain far below one nanosecond even 
for doubly-magic superheavy nuclei above A > 300. 
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FIG. 5: Potential strength parameter A (upper part) and 
volume integrals Jr (lower) versus decay energy E for su- 
perheavy nuclei. Known nuclei are shown as full circles. 
The extrapolated doubly-magic nucleus with Z mag i c = 132, 
A'magic = 194, and A mag i c = 326 is shown as open circle. The 
open squares are nuclei with Q = 18 and Q = 20 (see also 
Fig. . 
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FIG. 6: Potential strength parameter A (upper part) and 
volume integrals Jr (lower) versus decay energy E around the 
shell closure N = 82 for 13 °- 136 Xe <g> a (diamonds), 132 - 142 Ba 
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(squares) isotopes (see text; description of symbols see Fig.[4j|. 
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